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Abstract 

Background: Sensing and responding to ambient temperature is important for controlling growth and development 
of many organisms, in part by regulating mRNA levels. mRNA abundance can change with temperature, but it is 
unclear whether this results from changes in transcription or decay rates, and whether passive or active temperature 
regulation is involved. 

Results: Using a base analog labelling method, we directly measured the temperature coefficient, Qio, of mRNA 
synthesis and degradation rates of the Arabidopsis transcriptome. We show that for most genes, transcript levels are 
buffered against passive increases in transcription rates by balancing passive increases in the rate of decay. Strikingly, 
for temperature-responsive transcripts, increasing temperature raises transcript abundance primarily by promoting 
faster transcription relative to decay and not vice versa, suggesting a global transcriptional process exists that controls 
mRNA abundance by temperature. This is partly accounted for by gene body H2A.Z which is associated with low 
transcription rate Qio, but is also influenced by other marks and transcription factor activities. 

Conclusions: Our data show that less frequent chromatin states can produce temperature responses simply by virtue 
of their rarity and the difference between their thermal properties and those of the most common states, and 
underline the advantages of directly measuring transcription rate changes in dynamic systems, rather than inferring 
rates from changes in mRNA abundance. 



Background 

The mechanism for ambient temperature sensing in 
plants is unclear. Control of transcript levels is believed 
to be important in responses to temperature [1-4] but 
affects of ambient temperature on transcription and 
mRNA decay rates have not been measured. According 
to the work of Arrhenius [5] the temperature coefficient 
(Qio) of biochemical reactions is expected to be 2 to 3 at 
biological temperatures: yet less than 2% of Arabidopsis 
thaliana genes have a two-fold or greater difference 
in expression level between 17°C and 27°C [6]. The 
remaining genes either have rates buffered against chan- 
ging temperatures, or passive increases in transcription 
rate must be offset by a balanced increase in decay rate, 
leading to higher turnover but static steady state levels. 
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Despite this fundamental uncertainty, steady state tran- 
scriptomic responses to ambient temperature have been 
used to infer a role for chromatin modifications in 
temperature signaling [2,7]. 

4-Thiouracil (4SU) is a non-toxic base analogue that 
has been shown to be incorporated into mammalian and 
yeast mRNA during transcription [8-12]. Biotinylation 
and column separation allow 4SU-labeled RNA to be 
separated from unlabeled RNA, and transcriptomic ana- 
lysis using the separated samples can be used to simul- 
taneously calculate mRNA synthesis and decay rates [8]. 
Here we use 4SU labeling to measure transcription rates 
and determine the Qio genome-wide of mRNA synthesis 
and decay rates in Arabidopsis thaliana. We show 
that ambient temperature has large passive effects on 
both mRNA synthesis and decay rates, and that where 
temperature controls transcript abundance it does so by 
regulating transcription relative to decay and not vice 
versa. Our analysis suggests that transcription factor 
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binding sites and epigenetic state combine to create a 
complex network of temperature responses in plants. 

Results 

Cells incorporate 4SU into RNA and this has been 
exploited in mammalian cells [8,11,12] and in yeast [13] 
to measure mRNA synthesis and decay rates. In order to 
determine whether plants can take up 4SU we floated 
intact seedlings in MS medium and monitored 4SU in- 
corporation into RNA by biotinylation and dot blot (Fig- 
ure Sla in Additional file 1). This clearly showed that 
plants incorporate 4SU from the environment into RNA 
and that concentrations as low as 1 mM lead to a signal 
detectable above background within 1 hour (Figure IB). 
The resulting RNA could be separated from unlabeled 
RNA by biotinylation and passage through a streptavidin 
column as described previously. At 1.5 mM the flow- 
through can be depleted of detectable 4SU-labeled RNA, 
whilst labeled plant RNA is highly concentrated in the 
fraction recovered from the column [8,13] (Figure Sic in 
Additional file 1). To maximize recovery we chose a low 
concentration of 4SU at 1.5 mM [8] as high labeling fre- 
quencies are known to lead to binding of fewer more 
frequently labeled transcripts to the columns and reduce 
recovery. At this concentration Arabidopsis plants 
treated with 4SU showed the same growth and survival 
as control plants (Figure S2a in Additional file 1), sug- 
gesting 4SU has low toxicity in plants, as in other organ- 
isms. Therefore, 4SU dynamics in Arabidopsis seedlings 
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Figure 1 Schematic of experimental design and workflow. 

(A) Experimental design. Seedlings were floated on MS medium 12 
hours prior to the experiment, which started on addition of 4SU. 
RNA was harvested concurrent with 4SU addition, and then either 1 
or 2 hours later, depending on the temperature. (B) Experimental 
workflow for harvested samples. Samples at time 0 are extracted 
and hybridized directly, whereas samples at 1/2 hours are split. Total 
RNA is hybridized directly and compared to samples at time 0 to 
give change in steady state level during the labeling period. Total 
RNA was also split into two fractions by biotinylation and column 
separation, each of which is also hybridized to give abundance data 
for transcripts synthesized before and after 4SU addition. 



resemble those described for other experimental sys- 
tems. Preliminary experiments showed that RNA turn- 
over was faster at 27°C compared to 12°C (Figure S2b in 
Additional file 1), suggesting that temperature generally 
affected transcription rates. 

We designed an experiment to determine the Qio of 
mRNA synthesis and decay rates genome-wide. mRNA 
abundances were analyzed at two time-points at two 
temperatures 10°C apart, by microarray (Figure 1). At 
dawn total RNA was collected from seedlings floating on 
MS medium at 17°C and at 27°C, and 4SU was added to 
the remaining samples. After 1 hour at 27°C or 2 hours 
at 17°C material was flash-frozen, and RNA from this 
latter collection was separated into 4SU-labeled and un- 
labeled fractions (Figure Sic in Additional file 1). Total 
RNA from both the beginning and end of the labeling 
period, as well as labeled and unlabeled fractions at the 
end of the labeling period, were analyzed by Affymetrix 
AThl Genechips. Using linear regression [8] the relative 
contributions of the 4SU-labeled RNA (newly synthe- 
sized) fractions and unlabeled RNA (pre-existing, was 
transcribed before 4SU addition) to the total RNA frac- 
tion could be calculated, and the expression level of each 
gene normalized accordingly [8] (Figure S3 in Additional 
file 1; Materials and methods). Combining this informa- 
tion with the change (if any) in steady state levels of 
each transcript over time, from 12,997 genes called 
present by MASS we could calculate transcription and 
mRNA decay rates for 7,291 genes expressed in whole 
seedlings at both temperatures [8] (Materials and 
methods): dropouts included genes with highly variable 
expression (often low expressed) at one or both temper- 
atures, or genes with long half lives for which our label- 
ing period was too short to resolve. However, short 
labeling times are desirable in plants where transcrip- 
tional regulation can be highly temporally dynamic even 
in constant conditions. 

Plant transcription rates showed a log normal distribu- 
tion at both temperatures (Figure 2A; Additional file 2). 
The wide range of measured rates is good evidence that 
4SU incorporation into mRNA depended on variable tran- 
scription rates, rather than a single rate of uptake. Gene 
Ontology (GO) categories, including protein synthesis, 
metabolism, histones and transport, showed consistently 
high transcription rates compared to the whole genome, 
suggesting that genes with house-keeping functions have 
high transcription rates (Figure 2B,C). Categories with low 
rates included cell cycle gene expression, which presum- 
ably reflects the relative lack of cell division in mature 
plant tissues. The distribution of transcription rates shifted 
with temperature, showing that temperature has a major 
affect on absolute transcription rates in plants (Figure 2A). 

The calculated half-lives show wide variation between 
different individual genes (Figure 3A,C; Additional file 2). 
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Figure 2 Genome-wide distribution of transcription rates at 27°C and 17°C. (A) Histogram to show the distribution of transcription rates, in 
microarray units per hour, for 7,291 genes for which half-lives and transcription rates could be calculated. (B) Functional categories based on Gene 
Ontology terms of the 5% of genes with fastest and slowest transcription rates, at both temperatures. (C) Pie charts to show the distribution of 
transcription rates for all genes analyzed, at 17°C and 27°C compared to those encoding ribosomal proteins. Data points are calculated from the 
mean values for each locus from three biological replicate experiments. 
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Figure 3 Genome-wide distribution of mRNA decay rates at 27°C and 17°C. (A) Histogram to show the distribution of half-lives for 7,291 
genes for which half-lives and transcription rates could be calculated. (B) Functional categories, based on GO terms, of the 5% of genes with 
fastest and slowest decay rates (shortest and longest half-lives), at both temperatures. (C) Pie charts to show the distribution of half-lives for all 
7,291 genes, at 17°C and 27°C. (D) Pie charts to show the distributions of half-lives for genes encoding ribosomal proteins. (E) Pie charts to show 
the distributions of half-lives for genes annotated as encoding transcription factors. (F) Pie charts to show the distributions of half-lives for 
intron-less genes, at both temperatures. Data points are calculated from the mean values for each locus from three biological replicate determinations. 
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At 27°C the mean half Ufe for the genes analyzed is 
1.87 hours and this increases to 4.96 hours at 17°C. 
GO term-based analysis showed that gene categories 
enriched in genes with slow decay rates encode proteins 
with housekeeping functions, including basic metabol- 
ism and protein synthesis (Figure 3B,D), in agreement 
with previous analyses [14,15]. Genes with fast decay 
rates include those involved in signaling, particularly 
transcription factors, with 50% showing half-lives under 
1 hour at 27°C and under 3 hours at 17°C (Figure 3e). 
Our transcription and decay rate data could be con- 
firmed by RT-PCR using longer labeling time series, 
showing that rate determination by microarray pro- 
duces reliable data that can be independently corrobo- 
rated, despite only measuring values at two time-points 
(Figure S4 in Additional file 1). 

We found no correlation between decay or synthesis 
rates and mRNA length, UTR length, uracil content, or 
GC content (Figure S5 in Additional file 1). This is con- 
firmation that the labeling is sufficient and not biasing 
recovery of newly synthesized RNAs towards those with 
higher uracil content. However, we found that genes 
without introns have shorter half-lives (Figure 3F), cor- 
roborating similar evidence from a study using tran- 
scriptional inhibitors [15]. No relationship between 
rates of decay and transcription was observed and fold- 
change in transcript level during labeling correlated 
with decay rate but not with transcription rate (Figure S5 
in Additional file 1). 

A key question is whether changes in transcript abun- 
dance with changes in ambient temperature are related to 
passive temperature affects or represent the result of an 
active temperature signaling process. Average decay and 
transcription rates are both faster at 27°C than at 17°C. At 
27°C the mean transcription rate for all genes analyzed is 
236.1 microarray units per hour; at 17°C it is 73.0. Like- 
wise, the mean half-life at 27°C is 1.87 hours and at 17°C 
it is 4.96 hours: this demonstrates that both rates are over 
three times faster at 27°C compared to 17°C (Figures 2 A, 
3 A and 4). Thus, there is no general mechanism in plants 
that buffers mRNA abundance against temperature 
changes, and increasing temperature therefore increases 
mRNA turnover by affecting both synthesis and decay 
rates. We calculated the temperature coefficients (Qio) 
for all 7,292 genes for both transcription and decay rates 
(Figure 4 A, B). These were found to follow a log-normal 
distribution, indicating that the mean Qio is influenced by 
some genes with particularly large Qiq. The mean Qio at 
17°C for synthesis rates for all genes analyzed is 3.57 and 
the median 3.40. For decay rates the mean Qio at 17°C is 
3.29 and the median 2.80. Both means are higher than the 
Qio of two- to three-fold predicted by the Arrhenius equa- 
tion for a passive increase in rates, suggesting the presence 
of active response processes, as well as passive changes. 



Because splicing has also previously been shown to 
play a role in the control of gene expression by 
temperature [16], we examined the affect of intron num- 
ber on the distribution of synthesis and decay rate 
temperature coefficients (Figure S6 in Additional file 1). 
We found no evidence that intron presence or number 
could influence Qio of either transcript synthesis or 
decay on a global scale, suggesting that temperature af- 
fects a limited number of QioS in this way, or that 
temperature-regulated splicing does not frequently affect 
abundance of all detectable RNA species from a locus. 

Ambient temperature has been shown to modify the 
transcript levels of approximately 2% of the Arabidopsis 
genome [2,6]. Increasing transcript levels at higher tem- 
peratures could be achieved by a faster transcription rate 
or a slower decay rate. Our dataset includes 417 genes 
whose transcript abundances are more than 1.5-fold and 
significantly higher at the higher temperature, and 458 
genes with transcript abundances that are more than 1.5- 
fold and significantly lower at the higher temperature. 
These two sets of genes are clearly separated by Qio for 
transcription rate, but not decay rate QioS (Figure 4C). 
Therefore, a Qio of transcription that diverges from the 
mean is associated with mRNA species that change in 
abundance with temperature. Importantly, most tran- 
scripts that are more abundant at 27°C still show increases 
in decay rates with increasing temperature (Figure 4C), 
and increases in transcript abundance occur because tran- 
scription rate QioS exceed decay rate QioS. 

We divided our dataset into three groups: genes whose 
levels rise at least two-fold with increasing temperature 
(warm genes (WGs)), genes whose levels stay roughly 
constant (unchanged genes (UGs)) and genes whose 
levels rise at least two-fold with decreasing temperatures 
(cool genes (CGs)). We then calculated the Qio of tran- 
scription and decay for each group (Figure 4D, E). Tran- 
scription rate QioS of the WGs are significantly and 
remarkably higher than those of the UGs, which in turn 
are higher than those of the CGs (Figure 4E). In con- 
trast, only small changes in decay rate QioS were noted 
between the three groups and these were dwarfed by 
those observed for transcription rates. Together, our 
data shows that large changes in transcription rate are 
important in determining transcript abundance changes 
with ambient temperature, and that, on average, this is 
achieved by overcoming smaller, predominantly passive 
changes in decay rate. 

Because we can clearly show that regulation of tran- 
scription but not decay rates is central to the ambient 
temperature regulation of gene expression, we explored 
the possibility that key transcription factors or epigenetic 
states play a role in this process. For instance, previously 
it has been suggested that ARP6 and the SWRl complex 
control ambient temperature responses by functioning 
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Figure 4 Increasing temperature affects transcript abundance by active promotion of transcription rates. (A) Histogram to show the Qio 
of decay rates at 17°C for all 7,291 genes analyzed; the mean Qio for decay rate is 3.29 and the median 2.81. (B) Histogram to show the Qio of 
transcription rates at 17°C for all 7,291 genes analyzed; the mean Qio for transcription rate is 3.57 and the median 3.40. A number greater than 1 
indicates higher transcription/decay rates at 27°C than at 17°C. (C) Scatter plot showing the Qio of transcription rate plotted against Qio of decay 
rate for genes whose transcript abundances are at least 1.5-fold increased (red), 1.5 fold-decreased (blue) or unchanged (black) with increasing 
temperature. (D) Box plot to show the Qio of decay rates for genes with transcript abundances that are higher, lower or unchanged at 27°C 
compared to 17°C. Boxes show the median and the 5th, 25th, 75th and 95th per percentiles for each group of QioS. (E) Box plot to show the Qio 
of transcription rates for genes with transcript abundances that are higher, lower or unchanged at 27°C compared to 17°C. Boxes show the 
median and the 5th, 25th, 75th and 95th per percentiles for each group of QioS. 



in the replacement of his tone H2A with H2A.Z at lower 
temperatures [2]. The relationship between other epigen- 
etic states and temperature responsiveness has not been 
clearly addressed. In addition, signaling pathways con- 
trolled by temperature-regulated transcription factors 
have been shown to play a role in temperature signaling, 
including heat shock factors, basic helix-loop-helix 
(bHLH) transcription factors, FLOWERING LOCUS C 
(FLC), ABA response element (ABRE) -binding factors, 
and circadian clock genes [1,4,17-19]. Because our data 
showed that control of transcription rate is important in 
RNA abundance responses to ambient temperature, we 
measured the frequency of each feature in gene groups 
with QioS ranging from <1 to <8, based on published Ara- 
bidopsis seedling chromatin immunoprecipitation data at 
22°C [20,21], or the presence of an implicated transcrip- 
tion factor binding site up to 1 kb upstream of the tran- 
scription start site using The Arabidopsis Information 
Resource (TAIR) 10 annotation (Figure 5A). 

This analysis revealed some striking findings (Figure 5A). 
High levels of gene body H2A.Z [20] were more frequent 
in genes with QiqS lower than 2, with frequency rising 
from 8% of genes with average QiqS to 35% of genes 
with QioS lower than 1. This supports previous data sug- 
gesting that H2A.Z has a role in ambient temperature- 
regulation of transcription [2], but suggests that gene body 
H2A.Z rather than presence at the transcription start site 
is important for differential expression in response to 



temperature changes. In contrast, genes that are depleted 
in gene body H2A.Z showed little bias in Qiq. Analysis of 
the Qio of genes whose expression is altered in arp6 mu- 
tants at 12°C [2] found that genes with higher than aver- 
age transcription rate QiqS were unaffected in abundance 
in arp6 mutants at 12°C, whereas genes with QiqS below 
the mean were lower in abundance in arp6 mutants com- 
pared to the wild type. This further supports the hypoth- 
esis that between 12°C and 27°C H2A.Z/ARP6 activity is 
associated with low transcription rate Qiq. For other 
marks, H3K27me3 was frequently associated with genes 
with both high and low QiqS and was depleted in the gene 
bodies of passively regulated genes. In contrast, other 
more common marks, including DNA methylation, 
H3K4me3, H3K36me3 and H3k36me2, were preferentially 
present in genes with average QiqS. Given that the major- 
ity of genes have near-average QioS, the high frequency of 
these states combined with common properties deter- 
mines the mean temperature response of transcription 
(see Discussion for details). Together our data show that 
H2A.Z and H3K27me3 play a role in the generation of 
more thermostable transcription rates (that is, with a low 
Qio), but that H3K27me3 repression is also more fre- 
quently associated with greater than average QioS. 

Of the transcription factor binding sites tested, clear as- 
sociations between Qio and the presence of FLC binding 
sites and heat shock elements (HSEs) could be detected 
(Figure 6). As expected, HSEs were more frequent in the 
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Figure 5 Epigenetic control of transcription rate Qio in Arabidopsis. (A) The frequency distribution of chromatin marks among genes with 
differing temperature coefficients of transcription. Data (from Coleman-Derr et ol. [20] and Luo et ol. [21]) show the percentage of genes in each 
Qio category with each epigenetic character. Numbers of genes in each Qio group are 0 < Qio < 1, 41; 1 < Qio < 2, 527; 2 < Qio < 3, 2,039; 3 < 
Qio < 4, 2,447; 4 < Qio < 5, 1,389; 5 < Qio < 6, 526; 6 < Qio < 7, 174; 7 < Qio < 8, 82. (B) Pie charts showing the distribution of QioS of transcription 
rates for genes increased, decreased or unchanged in arp6 loss of function lines, using data from Kumar and Wigge [2]. (C) Box plot to show median, 
50% of data and 95% of data distributions for transcription rate QioS for ARP6 down-regulated genes (D), up-regulated genes (U) and genes with 
ARP6-independent expression (N). Asterisks show that the transcription rate QioS for genes with increased or decreased abundance in arp6 are 
significantly different from the genes with ARP6-independent expression (P = 0.0003 for genes decreased in orp6 and for increased P = 0.0001). 




temperature coefficient of transcription at 17 C 

Figure 6 Transcription factor binding site frequency in 1 Icb 
region upstream of the ATG of all genes in dataset classified 
by transcription rate Qio. Numbers of genes in each Qi 0 group 
are: 0 < Qio < 1, 41; 1 < Qio < 2, 527; 2 < Qio < 3, 2,039; 3 < Qio < 
4, 2,447; 4 < Qio < 5, 1,389; 5 < Qio < 6, 526; 6 < Qio < 7, 174; 7 < 
Qio < 8, 82. Circadian binding sites were determined using data 
from Covington et al. [19]; PLC binding site data are available from 
Deng et ol. [22]. 



promoters of genes with high transcription rate QiqS, 
showing that HSEs confer increased transcription in re- 
sponse to warmer temperatures, even throughout the am- 
bient range. In contrast, PLC binding sites were depleted 
in genes with average QiqS and more frequent in the pro- 
moters of genes with QiqS below 2 or above 5. This sup- 
ports an important role for PLC in configuring the 
ambient temperature response [17], and suggests that PLC 
can act as a repressor or activator of transcription in dif- 
ferent circumstances. Morning elements (MEs), evening 
elements (EEs), and G -boxes (bound by PIP subfamily 
transcription factors) showed a similar trend in that their 
frequency was marginally higher in the promoters of 
genes with lower QiqS, but then increased in the highest 
Qio category. Interestingly, ABREs were also present at 
higher frequency in the promoters of gene with the high- 
est QioS. Together this suggests that high transcription 
rate QiqS can be conferred by the combined action of 
multiple transcription factors. 

Discussion 

Measurements of Qio of specific cellular processes 
are not common: usually these are limited to whole 
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physiology, such as growth or gas exchange [23]. Here 
we show we can directly measure transcript synthesis 
and decay rates genome-wide in plants and use these 
data to calculate the temperature coefficients of synthe- 
sis and decay, which have means of 3.6 and 3.3, respect- 
ively, at 17°C (Figure 3). Large but balanced passive 
temperature effects on transcription and decay rate are 
observed for the majority of genes, and these ensure that 
mRNA abundance is constant at different growth tem- 
peratures. This will necessarily result in increased 
mRNA turnover at warmer temperatures, which will 
have a large energy cost to plant cells [24]. Importantly, 
differential control of transcript levels by ambient 
temperature is imparted predominantly by temperature- 
dependent transcriptional effects: these are presumably 
active processes because they produce changes in overall 
levels by failing to match passive effects on decay rate by 
virtue of high or low QiqS. Thus, a mRNA whose level 
increases at warmer temperatures does so because the 
transcription rate Qio exceeds the passive increases in 
decay rate. 

In order to generate an ambient temperature signal at a 
locus it is therefore necessary that the thermal response of 
this locus must deviate from the mean response (Figure 4); 
however, it follows that such a response must be relatively 
rare, otherwise it will itself define the mean. Previously it 
has been suggested that epigenetic effects play a role in 
the temperature control of transcription rates, with a glo- 
bal role ascribed to H2A.Z and a role for H3K27me3 spe- 
cifically in the case of FLC regulation [2,25]. We found 
that the presence of less common repressive chromatin 
states, such as gene body H2A.Z at 22°C, is correlated 
with reduced transcription rate Qio between 17°C and 
27°C (Figure 5). In contrast, more common chromatin 
states, such as H3K4me3, H3k36me2/3 and cytosine 
methylation, are biased towards genes with average tran- 
scription rate QioS. This is consistent with a view that the 
thermal responses of common states define the mean re- 
sponse, by virtue of their high frequency and similar prop- 
erties. The opposite Qio associations of DNA methlyation 
and H2A.Z is consistent with their known antagonistic 
roles [26]. 

Our data do not support a role for ARP6-dependent 
processes in generating a global transcriptional response 
to warm ambient temperatures, because genes with high 
transcription rate QiqS are expressed at normal levels in 
arp6 mutants (Figure 5C). However, we found a clear as- 
sociation between high gene body H2A.Z at 22°C and 
low transcription rate Qio between 17°C and 27°C and 
this relationship is corroborated by the observation that 
genes with low transcription rate QiqS are specifically 
mis-regulated in arp6 mutants at 12°C. Because this as- 
sociation specifically involves genes down- regulated in 
arp6 mutants at 12°C [2] (Figure 5C), gene body H2A.Z 



reduces Qio by catalyzing transcription at cooler temper- 
atures, relative to warmer temperatures. Our data are 
consistent with a model in which chromatin states have 
distinct activation energies of transcription, and that 
combining them in specific abundance ratios can pro- 
duce an ambient temperature response, much like com- 
bining two metals into a bimetallic strip can be used 
to create a thermometer. Alternatively, a 10°C rise in 
temperature could cause a mass rearrangement of mul- 
tiple chromatin marks that precipitates Qio variation. In 
principle, either could be used to generate a biological 
thermometer, but evidence is required that this process 
is transduced by plants into a temperature signal, as has 
been proposed during flowering time control or seed de- 
velopment [3,7]. The chromatin state data used for our 
analysis were chosen because the analyses were carried 
out at 22°C, equidistant between 17°C and 27°C. If 
temperature affects the frequency or location of less 
common marks, then this may also add an additional 
layer of complexity. 

H3K27me3 presence at 22°C is associated preferen- 
tially with extreme QioS across this temperature range 
(Figure 5A). Further work is necessary to understand 
whether it has a general role in transducing temperature 
signals. In the case of FLC, H3K27me3 is known to be 
required to stabilize transcription states but may also 
have a role in signal transduction [25,27,28]. 

Conclusions 

Our work reveals that transcription and mRNA decay 
are subject to both passive and active rate changes by 
temperature in A. thaliana. Our analysis reveals an 
interaction between chromatin state and transcription 
rate Qio such that chromatin marks bias transcription 
rate Qio relative to the mean, in particular gene body 
H2A.Z, which reduces Qio, and H3K27me3, which is as- 
sociated with extreme temperature responses. These 
interact with transcription factor activity to confer the 
observed ambient temperature response of the Arabi- 
dopsis transcriptome. Our approach shows that mea- 
surements of real rates have clear advantages over 
inferring rate changes from steady state levels, and will 
be useful for investigating dynamic plant responses to 
other signals. 

Materials and methods 

Plant growth and labeling 

Arabidopsis plants, ecotype Columbia, were grown in 
cycles of 12 hours 80 (imol white light, 12 hours 
dark. Eleven-day-old seedlings were removed from MS 
(Murashige and Skoog, Melfords, Ipswich, Suffolk, UK) 
agar plates and allowed to acclimatize overnight to float- 
ing in liquid MS media. They were then labeled by add- 
ing 4SU (Carbosynth Compton, Berkshire, UK) to the 
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liquid MS media 10 minutes after dawn. The concentra- 
tions added were different for different labeling times 
and temperatures and optimal concentrations for each 
set of conditions were determined experimentally using 
dot blots (see below for method) visualized using 
streptavidin-HRP (Genscript Piscataway, NJ, USA) to 
check the levels of 4SU in the RNA. Plant samples were 
removed from the liquid medium and flash frozen in li- 
quid nitrogen. 

RNA extraction and separation 

RNA was extracted from the plant samples using the 
RNeasy Plant (QIAGEN Germantown, MD, USA). The 
extracted labeled RNA was then biotinylated and sepa- 
rated on streptavidin columns as previously described 
[8]. The separated pre-existing RNA was tested by dot- 
ting onto Zeta probe membranes (Biorad Hercules, CA, 
USA) and detecting any remaining biotinylated RNA 
using streptavidin-HRP (Genscript), ECL plus reagents 
(GE Life Sciences, Little Chalfont, Buckinghamshire, 
UK) and a light detection system. RNA samples were 
checked for integrity using Agilent Bioanalyzer pico or 
nanochips, according to the manufacturer s protocols. 

Microarray analysis 

The RNA samples were labeled and hybridized to Affy- 
metrix ATHl chips using the manufacturer's protocols: 
these included three biological replicates for each sample 
type (total at end of labeling, labeled newly transcribed, 
unlabeled pre-existing and total at start of labeling) for 
each temperature (17°C and 27°C). Arrays were proc- 
essed using Mas5 with detection calls. Genes declared 
absent in one or more of the 12 total RNA samples were 
not included in further analysis. Transcript levels were 
scored as being significantly up- or down-regulated with 
temperature using SAM [29] with a false discovery rate 
of 1%. 

Normalization of newly transcribed and pre-existing RNA 
by linear regression 

Because equal amounts of RNA for the newly synthe- 
sized (NS), pre-existing (PE) and total (T) samples were 
loaded on the chips, measured transcript levels in the PE 
and NS samples were scaled to be directly comparable 
to the transcript level measured in the T dataset. To 
achieve this the ratios of the T mRNA accounted for by 
the NS and PE RNA pools were calculated by linear re- 
gression as described [8] based on the fact that for each 
gene the sum of the expression value ratios PE/T and = 
NS/T must equal 1 (100% of the total). But because 
equal amounts of RNA were loaded on the chips it is ne- 
cessary to use correction factors to get back to this value 
of 1. 



Correction factors that enable the transcript level of 
each gene in the NS and PE samples to be scaled to be 
directly comparable to those in the T sample were calcu- 
lated assuming: 

Cp.(PE/T) + Cn.(NS/T) = 1, (1) 

where Cp is the correction factor for pre-existing RNA 
and Cn is the correction factor for newly transcribed 
RNA. This rearranges to give: 

PE/T = (l/Cp)-(NS/T).(Cn/Cp) (2) 

To calculate the correction factors, NS/T for each 
gene is plotted against PE/T for each gene, giving a scat- 
ter plot and linear regression (Figure S5 in Additional 
file 1). The correction factors were calculated from the 
linear regression equation where 1/Cp is the intercept 
and Cn/Cp is the slope. 

Calculation of half-lives 

Half-lives were calculated assuming that mRNA decay 
follows a first order exponential decay kinetic [30,31]. 
For each gene, we can use the equation: 

tl/2 = -tL . ln(2)/ln(P/Ao) (3) 

where tl/2 is the half-life in hours, tL is the labeling 
time in hours, P is the scaled amount of pre-existing 
transcript remaining at the end of the labeling time and 
Aq is the amount of transcript at the start of the labeling 
time. 

Calculation of transcription rate 

The amount of RNA transcribed for each gene during 
the labeling period is equal to the amount of RNA de- 
tected at the end of the labeling plus the amount of NS 
RNA to have decayed during the labeling period. To get 
the amount synthesized in one hour we can use the 
formula: 



l_(g(-X tL)) 

where N is the scaled amount of newly transcribed RNA 
remaining at the end of the labeling time, tL is again the 
labeling time, and the decay constant X is: 

Ln(2)/tl/2. 
Transcript feature analysis 

cDNA and UTR length, GC content, U content and in- 
tron number as well as the presence or absence of short, 
potentially regulatory, sequences were extracted from 
the latest sequence files on the TAIR website [32] using 
Perl scripts (available on request). 
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Real-time RT-PCR 

RNA was labeled, extracted, separated and checked as 
above. First-strand cDNA was synthesized with 500 ng 
(or 100 ng for newly transcribed samples) of total RNA 
in 20 (il reactions. Superscript II Reverse Transcriptase 
(Invitrogen Carlsbad, CA, USA) and oligo dT (Invitro- 
gen), according to the manufacturers protocol Water 
(180 \A) was added and 2.5 \A of this was used in the 
quantitative PGR step. Primers are listed in Table S2 in 
Additional file 1. Quantitative PGR was carried out on 
the same samples for a number of genes and these were 
normalized collectively in a similar fashion to the micro- 
arrays, with the levels for the un-separated samples nor- 
malized to both UBQIO and TUB2. 

Data availability 

Microarray data generated for this manuscript are avail- 
able at the Gene Expression Omnibus under accession 
number GSE53071. Further datasets used in this analysis 
were GSE39045 [20] and GSE28398 [21]. 

Additional files 
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